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^ ; ABSTRACT 

c/3 , We use the combined GALFORM semi-analytical model of galaxy formation and GRASIL 

^ I ' spectrophotometric code to investigate the properties of galaxies selected via their sub- 

millimeter (sub-mm) emission. The fiducial model we use has previously been shown 

■ to fit the properties of local ULIRGs, as well as the number counts of faint sub-mm 
^> ' galaxies. Here we test the model in more detail by comparing the SEDs and stel- 

, lar, dynamical, gas and halo masses of sub-mm galaxies against observational data. 

■ We precisely mimic the sub-mm and radio selection function of the observations and 
Q!^ \ show that the predicted far-infrared properties of model galaxies with 6*850 > 5 mJy 

and SiA > SOfiJy are in good agreement with observations. Although the dust emis- 
0^ ' sion model does not assume a single dust temperature, the far-infrared SEDs are well 

I described by single component modified black-body spectrum with characteristic tem- 

00 . perature 32 ± 5K, in good agreement with observations. We also find evidence that 

' the observations may have uncovered evolution in the far-infrared-radio relation in 

ULIRGs out to z 2. We show that the predicted redshift distribution of sub-mm 
galaxies provides a reasonable fit to the observational data with a median redshift 
z = 2.0. The radio-selected subset of sub-mm galaxies are predicted to make up ap- 
proximately 75% of the population and peak at z = 1.7, in reasonable agreement with 
the observed radio detected fraction and redshift distribution. However, the predicted 
i^-band and mid-infrared (3-8/im) flux densities of the sub-mm galaxies (and LBGs) 
are up to a factor 10 x fainter than observed. We show that including the stellar 
TP-AGB phase in the stellar population models does not make up for this deficit. 
This discrepancy may indicate that the stellar masses of the sub-mm galaxies in the 
model are too low: M^, ~ 10^° M©, while observations suggest more massive systems, 
> 10^^ Mq. However, if the predicted K- and 3-8/im extinctions in the model could 
be dramatically reduced, then this would reduce, but not eliminate, this discrepancy 
Finally we discuss the potential modifications to the models which may improve the 
fit to the observational data, as well as the new observational tests which will be made 
possible with the arrival of new facilities, such as SCUBA2. 

Key words: galaxies: evolution - galaxies: formation - galaxies: high-redshift - sub- 
millimeter 



1 INTRODUCTION 

The discovery a decade ago of a population of faint submm- 
selected galaxies (SMGs) revolutionised o ur view of the cos- 
mic star formation history of the Universe (|Smail et al.ll 19971 : 



iHughes et al.lll998l : [B"arger et aLlligQel ). These galaxies (orig- 
inally discovered with the SCUBA instrument at 850//m) 
appear to be high redshift gal axies with star-formation 
rates exceeding lOOOMQyr^^ (,Ivison et al.ll2000l ; ISmail et al.l 
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bOOj ). These Ultra- Luminous Infrared Galaxies (ULIRGs) 
peak around z ~ 2 IChapman et all 12001 [ioosi ) and show 
a thousand fold increase in their abundance between 2 = 
and z = 2. \i this far-infrared emission arises solely from 
star-formation with a standard ("solar neighborhood") IMF, 
then they could potentially dominate the star-formation ac- 
tivity in the early Universe, dwarfing the contribution of 
galaxies selected in the rest-frame ultraviolet. The appar- 
ent intensity of these starbursts, the resulting high metal- 
licity, along with their large dynamica l masses, high gas 
fractions and inferred strong clu stering /Blain et alj 2004bl: 
Greve et al1l2005l:lTacconi et al.ir2006. 2008; Swinba nk et alJ 



20041 . 12OO6I ') are all suggestive of a close link to the forma- 



tion phase of the most massive s| 


Dheroids and black holes 


llLillv et al.ll 19991: ISmail et al.ll2002 


. 12004 IWebb et al. 


2003 


Genzel et al."2003': 'Smail et al.l 200 


4: Alexander et al. 


2003 


2005; Swinbank ct al. 2006). 



Reproducing the sub-mm galaxy p opulation has been a 
major challen ge for theoretical models l|Granato et al.|[2000l : 
iBaugh et al.l 12005). In part this is because (quite reason- 
ably) many of the recipes and constraints used to develop 
the models are based on the formation and evolution of "nor- 
mal" galaxies rather than the extreme populations. Hence 
these models have had difficulty reproducing extremely lu- 
minous galaxies with sufficient cool dust at high redshift 
without over-predicting the abundance of bright galaxies in 
the local Universe. Initial attempts to match the basic prop- 
erties of sub-mm galaxies (whilst maintaining the match to 
the present day Tf-band luminosity function and IRAS 60fim 
luminosity fu nction) under-pred icted the 850/xm counts by 
a factor 30 x (|Baugh et al.l l2005!) despite ACDM producing 
enough baryons in massive halos at 2; ~ 2 to match observa- 
tions of gas masse s in sub-mm galaxies l|Genzel et al.|[20o3 : 
iGreve et al.ll2005l ). 

One solution to th is problem w as to alter the Initial 
Mass Function ( IMF): iBlain et all l|l999l ) suggested that 
a Salpeter IMF (|Salpeterl 1 19551 ) "with a low-mass thresh- 
old of 3 M0 was required in far-infrared luminous galax- 
ies around 2 ~ 2 since the implied star-formation rate 
for a Salpeter IMF integrated to the canonical value of 
0.07 Mq would over-predict the integrated stellar density at 
2 = 0. Invoking a top-heavy (or "flat") IMF in the semi- 
analytic framework has been shown to provide a much im- 
proved fit to the sub-mm galaxy counts and redshift distri- 
bution as well as the Lyman-break galaxy luminosity func- 
tion whilst still maintaining a goo d fit to the propert ies of 
the present day galaxy populations (|Baugh et al.ll2005l . here- 
after B05). In this model, a standard IMF is adopted in 
quiescently star-forming galaxies, whilst the star-formation 
induced by galaxy mergers produces stars with a flat IMF: 
dn/dln{m) = with a; = (co mpared to a Salpeter IMF 
which has x = 1.35; ISalpeteill95"5l ). With a larger proportion 
of high mass stars, the energy radiated in the ultra-violet per 
unit mass of stars produced is increased, thus increasing the 
amount of radiation to heat the dust. Moreover, the flat IMF 
produces a higher yield of metals from type II supernovae, 
thus increasing the dust content of the galaxy and boosting 
the luminosity in the sub-millimeter waveband. 

This assumption of a flat IMF in bursts is controversial 
and so it is essential to explore the predictions of this model 
in more detail. In particular, the model can be (i) tested 
against the growing observational multi-wavelength data on 



sub-mm galaxies to test whether the choice of parameters is 
suitable and (ii) provide useful constraints on observational 
data, in particular for understanding possible selection bi- 
ases. For example, since much of what is known about sub- 
mm galaxies is based on the radio-detected sub-sample, one 
outstanding issue is how much this radio-selection (which 
does not beneflt from the negative K-correction experienced 
in the sub-mm waveband) has affected the conclusions be- 
ing drawn about SMGs. Are the radio-undete cted fraction of 
SMG s at either signiflcantly higher redshift ([Younger et al.l 
[2003) or do they instead hav e far-infrared colours w hich rep- 
resent much colder systems (|Chapman et al.ll2005^ ? 

In §2 we describe the basic properties of the galaxy for- 
matio n model we employ, but refe r the reader to ICole et al.l 
ll2000t).lBenson et al.f l|200^ . l2003t ). rBaugh et all (|2005l ) and 
iLacev et al. I (|2008l ) for a detailed description. In §3 we 
test the predicted far-infrared properties of sub-mm-selected 
galaxies against the available observational data. In §4 we 
discuss the masses and evolution of model SMGs compared 
to those estimated from observational data and in §5 we give 
our conclusions and future prospects for constraining the 
properties of sub-mm galaxies both theoretically and obser- 
vationally. We use a flat, ACDM cosmology with Qm — 0.3, 
as = 0.93 and Ho = 100ft kms'^ Mpc"^ with h = 0.7. 



2 THE MODEL 

The galaxy formation mod el whic h we e mplo y is the one 
descri bed in detail in .Baugh et al.l (|2005l ) and iLacev et al.l 
20081 ). The mod el is based on the semi-analytic code of 



Cole et al.l ||2000|) w ith im portant revisions described in 



Benson et al. 



Jwot I2003I ) and is reviewed extensively in 
Baughl (|2006l ). In summary, the formation, assembly and 
evolution of galaxies are calculated using a background cos- 
mology in which the dark matter structure grows hierarchi- 
cally. The physical ingredients considered in the model in- 
clude: (i) The formation of dark matter halos through merg- 
ers and accretion of material, (ii) The collapse of baryons 
into the gravitational potential wells of dark matter halos. 
(iii) The radiative cooling of gas that is shock heated during 
infall into the dark halo, (iv) The formation of a rotationally 
supported disk of cold gas. (v) The formation of stars from 
the cold gas. (vi) The injection of energy into the interstellar 
medium, through supernova explosions or the accretion of 
material onto a super-massive black hole, (vii) The chemi- 
cal evolution of the interstellar medium, stars and the hot 
gas. (viii) The merger of galaxies following the merger of 
their host dark matter halos, due to dynamical friction, (ix) 
The formation of spheroids during mergers due to the re- 
arrangement of pre-existing stars (i.e. the disk and bulge 
of the progenitor galaxies) and the formation of stars in a 
burst, (x) The construction of a star- formation history and 
a composite stellar population for each galaxy. 

The spectral energy distribution (SED) from this com- 
posite stellar populatio n is calculated usi ng the spectropho- 
tometric model GRASIL (ISilva et al.lll998l ). GRASIL computes 
the emission from both the stars and dust in a galaxy, based 
on the star formation, metal enrichment histo ry and the sizes 
as pr edicted by the semi-analytical model l|Granato et al.l 
I2OOOI ). GRASIL includes radiative transfer through a two- 
phase dust medium, with a diffuse component and giant 
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molecular clouds, and a distribution of dust grain sizes. We 
note that GRASIL does not assume a single dust temperature 
for the galaxy. Stars are assumed to form inside the clouds 
and then gradually migrate out. The output from GRASIL is 
the galaxy SED from the far-UV to radio wavelengths. 

2.1 Radio Emission and the Far-Infrared—Radio 
Correlation 

In addition to the sub-mm/far-infrared emission from dust, 
the radio emi s sion i s also included in the model following 
iBressan et al] l|2002h . The radio emission is produced by 
(i) thermal bremsstrahlung from Hll regions and (ii) syn- 
chrotron radiation powered by acceleration of relativistic 
electrons in supernova remnants. In local galaxies, there 
is a strong correlation between non-thermal radio and far- 
infrared emission, holding over five decades of luminos- 
ity (e .g. iHelou et al.lll985l : lYun et al.ll200ll : IVlahakis etahl 
I2OO7I ). The standard explanation of this relationship is that 
both the far-infrared and the bulk of the radio emission are 
caused by high-mass 5Mq) stars. These stars both heat 
the dust (which then emits far-infrared emission) a nd at the 
end o f their lives explode as supernovae (see e.g. ICondonI 
Il992l and references therein), producing the relativistic elec- 
trons responsible for synchrotron radiation. The strength of 
this synchrotron emission depends on not only the number 
density of relativistic electrons, but also the strength of the 
magnetic field. The rate of energy input into the relativis- 
tic electrons depends on the assumed low-mass cut off for 
Type II supernova produ ction, which is usually taken as 
8M0 (jBressan et al.H2002l ). 

In the model galaxies, both the thermal and non- 
thermal contributions to the radio emission are calculated. 
The thermal component of the radio emission is estimated 
from the stellar emission essentially without free parame- 
ters, assuming that all ionizing photons are absorbed by 
gas within the galaxy, and is proportional to the instan- 
taneous ionizing luminosity of the stars. The synchrotron 
(non-thermal) radio emission is assumed to be proportional 
to the instantaneous Type II supernova rate, but its calcu- 
lation involves two empirical parameters, which are essen- 
tially the efficiency with which the supernova explosion en- 
ergy is converted into energy of relativistic electrons, and the 
power-law index of the injection energy spectrum of these 
relativistic electrons, which determines the spectral index of 
the synchrotron emission. It is assumed that the radiative 
lifetime of the relativistic electrons is short enough that the 
total synchrotron luminosity is always equal t o the instanta- 
neous energy injection rate from supernovae. iBressan et al.l 
l|2002t l estimate both of these parameters from observa- 
tions of local galaxies; in particular, the efficiency factor 
was chosen based on observations of the Milky Way galaxy. 
They found that this choice also reproduced the observed 
ratio of radio to far-infrared emission for normal spirals. 
In this model, the ratios between the far-infrared and ra- 
dio luminosities and the instantaneous star formation rate 
(SFR) all approach constant values if the SFR varies only on 
timescales ^100 Myr (as in normal spirals). However, due to 
the different lifetimes of the stars powering the far-infrared, 
thermal and non-thermal radio emission, there are timelags 
between the corresponding luminosities and the formation 
of the stars involved, and so these ratios vary if the SFR 



changes on timescales ^lOO Myr, as is the case in s tarbursts. 
This is discussed in detail in IBressan et all (|2002l ). and has 
consequences for the results in this paper. 

In Fig. [T] we show the far- infrared-radio (60/im versus 
1.4 GHz) correlation for local star- forming galaxies and over- 
lay the model predictions for galaxies at z < 0.1. This shows 
that the model reproduces the form of the observed far- 
infrared-radio correlation at lower luminosities. However, at 
higher luminosities the p redicted radio-far - infrared relation 
lies below that obser ved (jHelou et al.ll 19851 : IYuu et al.ll200ll : 
Vlahakis et "all 120071 '). by a factor 3.2 ± 0.2x, assuming the 
Bressan et al.l (|2002 ) normalisation of the radio emission us- 
ing the Milky Way. As we are interested in the evolution of 
SMGs, which have ULIRG-like luminosities [Ltoi > lO^^L©) 
and appear likewise to be dusty starbursts, it seems appro- 
priate to renormalize the non-thermal radio emission in the 
model so that it matches the properties of high-luminosity 
IRAS-selected samples at z = 0. This seems likely to pro- 
duce the most realistic estimates for the radio luminosi- 
ties of SMGs, in the framework of the present galaxy for- 
mation model. Physically, this renormalization corresponds 
to assuming a higher efficiency for conversion of supernova 
blastwave energy into energy of relativistic electrons in star- 
bursts as compared to normal spirals, which might result 
from the higher densities or magnetic fields in starbursts. 
In the present paper, we have therefore multiplied the non- 
thermal radio luminosities from GRASIL by a factor 3.2, as 
shown in Fig. [T] We note that this renormalization might 
not be needed at 2: = in a galaxy formation model which 
predicts significantly d i fferen t burst timescales in ULIRGs. 
Indeed, IBressan et al.l l|2002l ) reproduced the SED of the 
ULIRG Arp220, including its radio emission, using their 
standard synchrotron normalization, but treating the burst 
timescale as a free parameter. However, in our galaxy forma- 
tion model, the burst timescales are already fixed by other 
considerations. 

In what follows, we parameterize the relative strength 
of the far- infrared and radio emission using ql, defined as 
qL = ioy(L6oi/([4.52THz] Li.4GHz)) where Lboi is the total 
dust luminosity and Li.4ghz is the rest-frame radio lumi- 
nosity. The latter is calculated from the observed 1.4GHz 
fiux 51.4 using 1/1.4 GHz = 47rZ)|,Si.4(l + z)°'~^, where the 
spectral index, a, has a typical value of 0.7 for non-thermal 
sources; ICondonllT993 ). For local ULIRGs, ql is measured 
to be qL = 2.34 ± 0.01 with a^^ = 0.10 (e.g. lYun et all 
[200I). As described above, we have normalised the radio 
emission of 2 < 0.1 model ULIRGs so that by construc- 
tion they have qL = 2.34. The model predicts a disper- 
sion — 0.11. We stress that this renormalisation en- 
sures that our model agrees with observations of the radio 
and far-infrared emission of galaxies in the luminosity range 
appropriate for SMGs in the local Universe. However, the 
model also predicts evolution in these quantities. For an ex- 
ponentially declining star-formation rate with a top-heavy 
IMF, the radio ffux density will plateau shortly after the ini- 
tial burst since the SNe are produced from stars with mass 
^ 8 Mq which have lifetimes of approximately 50 Myrs. In 
contrast, the far-infrared luminosity will continue to increase 
for longer as dust is produced and then heated by new gen- 
erations of stars extending down to lower masses ^ 5 Mq . 
As a result, the value of qL is not constant for bursts, but 
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Figure 1. Far-Infrared — radio correlation for galaxies in the 
GRASIL model compared to the SCUBA Local Ga l axy S urvey 
(SLUGS) and IRAS 2 Jy survey from lYun et alj 1 I2OOII) and 
IVlahakis et al ] l|2007l l. The red error-bars show the scatter in 
the observed data. The shaded area denotes the approximate 
60Atm and ra d io lu minosities encompassing the ULIRGs from 
ISanders et al. 1 1I2OO3I '). The model {z < 0.1) predictions are over- 
laid for both the original radio luminosities (open circles) and 
after the renormalisation factor 3.2x has been applied to the ra- 
dio flux densities (filled circles). We adopt this renormalisation 
of the model to ensure that it reproduces the far-infrared and 
radio properties of local galaxies with luminosities comparable to 
those observed for high-redshift SMGs — giving us a reliable local 
benchmark. 



depends both on the S FR timescale and on the age at which 
the burst is observed (jBressan et al.|[2002l ). 

The average qL in bursts therefore depends both on 
redshift and on how the galaxy sample is selected. For model 
ULIRGs a.t z = 2 (which overlap substantially with SMGs), 
we find qL ~ 2.19 ± 0.05 and Uar = 0.26 . Moreover, the e- 



of sub-mm galaxies comes from (jChapman et al.ll2005l here- 
after COS) who surveyed a total area of approximately 0.25 
degree^ across seven fields, securing spectroscopic redshifts 
for 73 sub-mm galaxies with 850/im flux densities > 5pJy 
and radio counterparts with flux densities 5*1.4 > 30/iJy. 
This comprises the primary comparison sample for our anal- 
ysis. 

The C05 survey relies on radio identification as the only 
secure method to pin-point the counterparts of large sam- 
ples of the sub-mm emitting galaxies. This is because the 
large beam-size (~ 15" for the JCMT at 850/im) conspires 
with the large number of possible counterparts making it 
impossible to identify the galaxy responsible for the far- 
infra r ed emission with 8 50/im imaging a l one: llSmail et ahl 
l200d : llvison et al.lll998l . I2OO0I . l2002i . |2005| . boOTh . Using the 
deepest radio imaging currently available (cr ~ 5-10/iJy), 
this radio selection identi fies ^ 65-80% of the brigh t (Sgso > 
5 m Jy) sub-mm galaxies (jlvison et al.l |2005| . |2007| ) , but po- 
tentially introduces biases when extrapolating the properties 
of radio identified sub-mm galaxies to the whole population. 
We investigate the effect of this bias in the model by defin- 
ing two classes of sub-mm galaxy: first, all sub-mm galaxies 
with Sgso > 5mJy are called SMGs; second, we define the 
subset of this population which are detectable in the radio 
(with Sgso > 5mJy and Si. 4 > SO^Jy) as radio-identified 
sub-mm galaxies or rSMGs. 

In order to simulate observational sub-mm catalogues, 
we must model the measurement noise. In particular, since 
most sub-mm catalogues are cut at a threshold around 3.5- 
4(7, sub-mm maps sufTer from flux boosting in which the flux 
limit of low significance sources can be increased above the 
survey signal-to-noise due to (a) confusion (the contribution 
of fainter sources within the large beams) and (b) the inclu- 
sion of low-significance sources which are boosted above the 
survey signal-to-noise due to the coincidence alignment with 
positive noise spikes. In order to properly compare the ob- 
servations with the models we therefore convolve the model 
sub-mm and radio fl ux densities with a typical 1-a noise o f 



folding timescale for star formation in model ULIRGs more o"S50 = 1.5 mjy (e.g. IScott et al. ] |2002l : ICoppin et al 



than halves between z = Q and z = 2 {Te/f = 300_i5g Myr 
at z — compared to 130lgQ° Myr at 2 = 2). In addition, the 
median age of the burst in model ULIRGs at z = 2 is much 
shorter than at z = (8lg° compared to 5014°° Myr). The 
young ages of the high-redshift model ULIRGs, combined 
with the time-lag between the onset of the burst and the 
production of the first supernova, thus results in an evolution 
in the far-infrared-radio relation for these luminous galaxies. 
We will return to this point when studying the radio and 
far-infrared properties of SMGs later. 



3 ANALYSIS 

3.1 Simulating Sub-mm Catalogues 

To test the GALFORM / GRASIL model we need to compare its 
predictions to the detailed multi-wavelength properties of 
high-redshift sub-mm galaxies. To achieve this we have to fo- 
cus on samples with reliable galaxy identiflcation (which has 
so far relied on high resolution radio imaging), and follow- 
up spectroscopy to provide precise and unambiguous red- 
shifts. By far the largest and most secure redshift survey 



Ivison et al 



20071) 
20021 . 



and 0-1.4 = 5-10/iJy respectively (e.g 

|2007| ). thus allowing a like-for-like comparison with obser- 
vations. This measurement noise is assumed in all following 
sections and analysis, and is included in all of the estimates 
which we give for average masses and luminosities of model 
galaxies. 



3.2 Number Counts, Redshift Distribution and 
Space Density of SMGs 

In Fig. [2] we show the predicted cumulative source counts 
as a function of 850/im flux density compared to the ob- 
servations derived from a number of different surveys with 
SCUBA. The fiducial model of B05 agrees well with the 
observations down to the deepest fiux density limits, with 
the model counts above 0.1 mJy being dominated by ongo- 
ing merger-driven bursts at high-redshift, as these authors 
demonstrated. 

To determine the effect of the radio-preselection on the 
number counts, we also show the cumulative source counts 
for galaxies which also have radio flux densities 5*1.4 > 
30^Jy. At a characteristic 850/im flux density limit of 5 mJy 
of C05, the model suggests that the radio detected fraction 
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Figure 2. The cumulative number counts per degree^ at 850/im 
as a function of sub-mm flux density. The model predictions are 
shown for SMGs and rSMGs respectively. We also plot the obser- 
vational re sults fro m various surveys with SCUBA on th e JCMT 
llSmail et al., ,2002; Cowie et al.ll2002l : ICoppin et ai1l2006l) . Above 
5mjy, the model suggests that the surface density of sub-mm 
galaxies with radio flux densities above 30/iJy is 75%. This agrees 
well with the observed radio-detected fracti on of Ssso >5mjy , 
Si. 4 >30/.tJy sub-mm galaxies of 65-80% l llvison et al.l |2005|. 
l2007t) . 
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Figure 3. The predicted redshift distribution of model SMGs 
and rSMGs. Left: N{z) for galaxies selected at both sub-mm and 
radio- wavelengths. The median redshift in the model rSMGs is 
z = 1.7 with (Tz = 0.8. We o v er-plo t the redshift distribution of 
rSMGs from [Chapman et al.l ||2005| ) (corrected for spectroscopic 
incompleteness in the redshift desert). Right: N{z) for a pure sub- 
mm selected sample of galaxies with 850/im flux densities greater 
than 5 mjy, peaking at z = 2.0 with riz = 1.0. Again w e over-plot 
the redshift distribution from IChapman et al.l (|2005| ) , corrected 
for spectroscopic incompleteness and radio incompleteness using 
a simple model for the evolution of the radio luminosity function. 
In both panels, the redshift distribution for the model galaxies is 
slightly shallower than the observations suggest (by Az ~ 0.2- 
0.3), although since the field-to-field variance in the observations 
is Az ~ 0.25 we conclude that (within the observational uncer- 
tainties) the redshift distributions arc in good agreement. 



is ~ 75%, consistent with the fraction typically f ound in 
sub-mm surveys (65-80%) (|lvison et al.ll2005l . |2007| ). 

The redshift distribution of the radio-identified SMGs 
from C05 (crudely corrected for spectroscopic incomplete- 
ness due to the redshift desert) is well fit by a Gaussian 
profile with a median redshift z = 2.0 and CTz = 0.7. Using a 
simple model to account for radio incompleteness, C05 esti- 
mate the underlying redshift distribution of SMGs to peak 
at z = 2.2 with a 1-a width of az = 1.3. 

In Fig. [3] we show the redshift distribution of SMGs and 
rSMGs in the fiducial model of B05 compared to observa- 
tions. The model SMGs have a median redshift of < z 2.0 
with az — 1.0. However, as Fig. |3] also shows, the radio- 
identified sub-set peaks at < z 1.7 with az — 0.8. Al- 
though there is reasonable overlap between the model and 
observations, the model SMGs and rSMGs both appear to 
peak at slightly lower redshift than inferred observationally 
(Az ~ 0.2-0.3 in both cases). We note that the field-to-field 
variation between the seven sub-fields in the C05 sample is 
Az ~ 0.25. This may be a si mple refiection of t he fact that 
SMGs are highly clustered (|Blain et al.ll2004bl ). and so we 
expect that the redshift distribution from C05 is likely to be 
uncertain by at least this amount due to cosmic variance. 
Therefore we conclude that the model and observations are 
in good agreement. 

It is also worth noting that IClemente et all (|2008l ) de- 
rive a median redshift z ~ 1.5 using photometric redshifts 
(rising to z ~ 1.9 when including the far-infrared photom- 
etry) for a similar sample of radio and mid-infrared identi- 
fied SMGs from the S HADES submm survey. However, as 
[elements et al.l (|2008l ) note, there are large errors on pho- 
tometric redshifts for sub-mm galaxies (Az 0.5-1), espe- 
cially for the sources with the faintest counterparts. Nev- 



ertheless, their results are also consistent with the model 
redshift distributions for the rSMGs. 

Using the number counts and redshift distribution, we 
can also compare the space densities. The observed space 
density of rSMGs from C05 indicates that between z = 0.9- 
3.5 the volume density should be 8.0 x 10"^ Mpc"^. In 
comparison, the predicted space density for model rSMGs in 
the same redshift interval is 1.1 ±0.1 x 10~^ Mpc~^, slightly 
higher but consistent with the observations. In Fig. [4] we 
show the space densities over the redshift interval z = 0.9- 
3.5 for the model SMGs compared to the SMGs in the C05 
sample. To interpret this data as a cumulative luminosity 
function, we also convert the number counts and 850/.im flux 
densities to approximate space densities and bolometric lu- 
minosities assuming a median redshift z = 2.0. 



3.3 Far-infrared Spectral Properties 

Constraining the far-infrared SED of sub-mm galaxies is a 
key observational goal since the thermal emission from cold 
dust dominates the bolometric luminosity of these galax- 
ies (Lg-iooOfim ^ 0.95L(,oi)- Sub-mm photometry (at wave- 
lengths other than 850pim) can be used to constrain the 
SED, measuring the apparent temperature of the dust, and 
infer their far-infrared luminosity and thus star-formation 
rate (assuming an IMF). These diagnostics also allow these 
galaxies to be placed in the context of other populations of 
high-redshift star-forming galaxies and active galactic nuclei 
(AGN). 

To test how well the SEDs predicted by the model galax- 
ies reproduce the dust emission spectra of SMGs, we com- 
pare the predicted far-infrared colours with existing obser- 
vational data. 
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Figure 4. The cumulative number counts per degree^ at 850/im 
for galaxies between z = 0.9—3.5. The model predictions including 
all redshifts are shown by the solid line. We also show the expected 
counts for radio-flux density cuts of Si. 4=0 and 30/iJy, respec- 
tively, only considering galaxies in the range z = 0.9—3.5. These 
are compared to the observational counts for SMGs from the C05 
redshift survey which are radio selected with Si. 4 ^ 30^Jy, and 
also in the range z = 0.9-3.5. The right-hand axis denotes the 
approximate space density in physical units, whilst the top axis 
denotes the approximate bolometric luminosity for galaxies at 
2 = 2.0 given their sub-mm flux density. 



3.3.1 1200/850/1171 Colours 

First, we compare the 850/im flux density ratio with 
longer wavelength photometry. In particular, it has been 
claimed that the millimeter /sub-millimet er flux density ra- 
tio is sensit i ve to redshift above z = 3 (|Eales et al ] |2003l : 
iGreve et all l2004l : llvison et all l2005h . We therefore com- 
pare the far-infrared photometry predicted by the model 
SMG s and rSMGs wit h recent observational constraints 
from iGreve et al.l (|2008l ) who surveyed the GOODS-North 
region at 1200/im with the Max-Planck bolometer array 
(MAMBO). For an 850/im selected sample, the model 
galaxies have a median 850/1200/im flux density ratio of 
5860/51200 = 2.4±0.2 between z = 1.5-3, and 5850/51200 = 
2.2 ±0.2 for a 1200/im selected sample with (T1200 = 0.7 mjy 
over the same redshift range (we note that the error-bars 
denote the 17 and 84%-ile range, which corresponds to the 
1-0" scatter for a Gaussian distribution). Both the 850- and 
1200/im-selected samples show very little evolution in colour 
over the redshift range z — 0.5-5 and are in good agreement 
with observational constraints (Fig. [5]). 



3.3.2 350/850/im Colours 

A better test of the far-infrared properties comes from 
shorter wavelengths (at 2 = 2 a black-body with a charac- 
teristic temperature of Td=35K peaks at an observed wave- 
length of ~ 300/im). Prior to the launch of Herschel, the 
most promising route to constrain shorter wavelengths in 
the far-infrared are 350- and 450/im photometry where the 
atmospheric transmission allows the brightest SMGs to be 
detected in good conditions. 

In Fig. Owe compare the predicted Ssso/Ssso colours 



with observational constraints. At 2 ~ 2 the median flux 
density ratio for model SMGs is 5350 / 5*850 = 2.5 ± 2 . 0. For 
comparison, iKovacs et al.l (|2006l ) and ICoppin et al.l (|2008| ) 
study a total of 27 SMGs at 350/tm which have secure 
redshifts. These SMGs have a median flux density ratio 
535o/5850 = 4.012 2 (this value includes upper limits for 
>^ 350/tm flux densities on three non-detections). Since the 
^ data at 350 and 450/im suffer from both low number statis- 
stics and low signal-to-noise detections we bin the data into 
f two redshift bins, z < 2 and z > 2. As this figure shows, the 
model predictions seem broadly consistent with the obser- 
vational data. There is a hint (< 1.5ct) that the model pho- 
tometry and data have opposite trends with redshift, but 
firm conclusions cannot be drawn with the current observa- 
tional data. However, we stress that this potentially power- 
ful test will be significantly improved with observations of 
larger samples at higher signal-to-noise from upcoming sur- 
veys with SCUBA2 and Herschel. In the meantime we note 
that the model predicts an evolution in the Saso/Ssso flux 
density ratio for the SMGs and rSMGs is well described by 
S350/S850 = -0.6 + 11.6 X (1 + 2)-i-^^ 

A more useful comparison of the observational and 
model SEDs can be made by comparing the bolometric lu- 
minosities estimated from fits to both real and model pho- 
tometry. To do this, we fit the 350- and 850/im photome- 
try with a modified blackbody of the form oc Bv{Td)i'^ 
where Bu is the Planck f unction evaluated a t the emitted 
frequ ency v, and /3 = 1.5 (jPunne et al . 2 00^ : ICoppin et al] 
l2008h . The result of this fit is a characteristic dust emis- 
sion temperature Td for each galaxy, and also an estimate 
of its bolometric dust luminosity Lb oi. Kovacs et al.l (|2006l ) 
(see also more recentlv lCoppin et al.ll2008l ) demonstrate that 
sub-mm galaxies have characteristic temperatures consistent 
with local starbursts (T ~ 35 K) , but are at least an order of 
magnitude more luminous in L(,o!. We can apply the same 
fitting procedure to the 350- and 850/tm photometry for all 
our model SMGs. In this way, we find a median characteris- 
tic emission temperature for both model SMGs and rSMGs 
of Td — 32 ± 5K. (Note that simply a characteristic tem- 
perature, since GRASIL does not assume a single dust tem- 
perature even within a single galaxy.) We also find that the 
single modified black-body fit (with /3 = 1.5) yields an esti- 
mated bolometric luminosity which on average is 0.94±0.25 
of the true model value (integrated between rest-frame 8 and 
1000/im) . We find that the median bolometric luminosity for 
model SMGs is 2.0 ± 1.5 x 10^^ Mq. Hence it appears that 
the far-infrared SEDs of model SMGs are broadly similar to 
what is inferred from the observations. 

With the bolometric luminosities of the model SMGs 
in hand, we can also investigate how the star-formation 
rates compare to those inferred observationally. The main 
difference in the star-formation rates comes from the 
adoption of the fiat IMF in bursts. This results in a 
much lower star-formation rate per unit bolometric lu- 
minosity for SMGs (which are dominated by merger in- 
duced bursts). We find that the relation between bolo- 
metric luminosity and instantaneous star-formation rate is 
SFR(O.15-125M0)(M0yr-i) = 1.01 x 10^'''' (ergs^^)) 
which means that an average model SMG with L(,oi = 
2 x 10^^ L© has a star-formation rate of SFR=77M0yr"^ 
compared to a typically observationally-derived value of 
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SFR(0.1-100 Mb )= 4.5x IQ-^^ Lg-iooo^m (ergs-^) using the 
iKennicutd calibration for a Salpeter IMF. 



3.4 Radio properties 

Although the 1.4 GHz radio emission is used to identify 
the galaxy responsible for the sub-mm emission, the sub- 
mm/radio flux density ratio has also b een used as a diagnos - 
tic of both temperature and redshift (|Carilli fc Yun|[l999l ). 
However, since the far-infrared colours scale with z) /Td, 
redshifting a fixed SED template has the same effect as 
changing the temperature at a fixed redshift, and so without 
knowledge of the temperature, far-infrared colours cannot be 
unambiguously used to derive redshifts. Indeed, with a se- 
cure redshift for a sub-mm galaxy, the effects of radio identi- 
fication become apparent: a canonical 5mJy radio-identified 
SMG at z = 2.4 with a 50/xJy radio counterpart has a char- 
acteristic dust temperature of 32 K. An increase or decrease 
in the dust temperature of just 10 K has a dramatic (fac- 
tor 10 x) effect on the sub-mm and radio-fiux densities and 
therefore target selection: as the temperature increases for 
a fixed luminosity, the hotter SEDs mean that the 850/im 
flux density falls below the detection threshold of 5mJy. 
Similarly, the coolest SMGs at 2 = 2.4 would have lower 
radi o-flux densities making an SMG undetectable in the ra- 
dio l|Chapman et al. | |2004 -1 iBlain et al. I l2004bl . l2003 ). 

In Fig. [5] we compare the predicted 5'85o/'S'i.4 flux 
density ratio. Both the observations and model galaxies 
show strong evolution with redshift (a result of strong K- 
corrections and evolution). The model SMGs have a me- 
dian flux density ratio S'85o/<S'i.4 = ISOlgg'' at z ~ 2-3, 
and so a model SMG with Sgso = 5mJy has 5*1.4 = 30/iJy. 
In the model, the S850/S1.4 for model SMGs evolves as 
S85o/Si.4 = -200 + 196 x(l + 

Whilst the model SMGs appear to reproduce the gen- 
eral evolutionary trend shown by the observations, over the 
redshift range z = 0.5-3.5 they also have a S85o/S'i.4 flux 
density ratio which is systematically too high compared to 
the observations, 1.26±0.24x (where the error is a bootstrap 
estimate), despite matching the 2 < 0.1 far- infrared-radio 
correlation for ULIRGs. For the model rSMGs, we derive 



the radio and submm flux sele ction do not i nfluen ce this 
value. In terms of observations, iKovacs et all (|2006l ) exam- 
ine the far-infrared-radio correlation for fifteen high redshift 
SMGs at 2 ~ 2 and derive < qL >= 2.14 ± 0.07 with an 
intrinsic spread cjq ~ 0.12. Compared to ultraluminous far- 



infrared galaxies at 2 = 0, which have qL = 2.34_ 



(see 



>j2.1|l . we see that the model predicts a similar level of evo- 
lution to that seen in the observations. More importantly, 
the fact that qL does not change for the model SMGs when 
we include the obser vational biases sugg ests that the low 
value of gL found bv lKovacs et al.l (|2006l ) is not a result of 
sampl e selection. Instead, it appears that the lKovacs et all 
(|2006 r) result reflects real evolution in the far-infrared-radio 
correlation. Whether this results from the same starburst 
age-related evolution found in the model, or from other pro- 
cesses, such as modest contributions to the observed radio 
flux densities from AGN activity, requires further observa- 
tional work. 
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Figure 5. Predicted sub-mm and radio flux density ratios for 
sub-mm galaxies as a function of redshift. Top: S850/S1.4 colours 
of sub-mm galaxies in the C05 sample compared to model SMGs 
and rSMGs. We plot individual galaxies in the C05 sample and 
show the median values in redshift intervals of = 0.5. We also 
plot the predictions for model SMGs and rSMGs. The error-bars 
indicate the the typical 67%-ile of the model predictions, and, for 
clarity are shown on alternate points. The model rSMGs have a 
systematically larger (~ 1.26it0.24x) S850/S1.4 flux density ratio 
than observations suggest over the redshift range z = 1-4. Mid- 
dle: S350/S850 colours as a function of redshift for model SMGs 
and rSMGs. The comparis o n sam ple i s taken from observat ions at 
350Atm from lKovacs et al.l l|2006l ') and lCoppin et al.l (|2008h which 
have good (> 3cr) detections at 350^m. Within the observational 
errors the model SMGs are in good agreement with observations 
(which suffer from both low number statistics and low-signal to 
noise detections). To show the redshift trend in the observational 
data more clearly, we bin the data into two redshift bins {z < 2 
and z > 2) (solid squares). Bottom: S850/S1200 colours as a func- 
tion of redshift for model SMGs and rS MGs compared to obser- 
vations of the GOODS-North regions bv lGreve et al.l l|2008l '). The 
solid lines denote the evolution for an 850/im selected sample, 
whilst the dotted line represents the evolution for a 2.5 mjy se- 
lected sample at 1200/im with (71200 = 0.7 mJy. Both the model 
SMGs, rSMGs and 1200/im selected model SMGs show remark- 



8 Swinbank et al. 



3.5 Restframe UV Properties 



3.7 Restframe Near-infrared Properties 



With precise positions and redsliifts available for the radio- 
identified SMGs in COS, this sample has also been an impor- 
tant basis for more detailed study of the morphologies, stel- 
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a combination of high resolution imaging with HST, as well 
as multi-band photometry and spectroscopy from optical, 
near- and mid-infrared wavelengths. In particular, the stel- 
lar populations have been probed using optical, near- and 
mid-infrared colours to constrain the stellar SED from the 
rest-frame UV to near-infrared. 

The next step in our comparison is therefore to exam- 
ine the stellar populations and masses of model SMGs. We 
start by comparing the predicted B- and ii-band magni- 
tudes of SMGs to those from C05. At the median redshift of 
the C05 sample, the observed B- and i?-bands sample the 
rest-frame UV, which is very sensitive to the instantaneous 
SFR, as well as the level of dust extinction. For model SMGs 
in the redshift range z = 0.9-3.5 the mean B- and _R-band 
magnitudes are Bab = 24.0 ± 1.3 and Rab = 22,.9t\ °o- In 
comparison the SMGs from C05 have Bab = 24.8 ± 1.2, 
Rab = 24.3 ±2.0, slightly fainter and redder than the mod- 
els predict but with significant overlap (just considering the 
median {B — R) colours, we note that an additional ob- 
scuration of E{B — V) = 0.2 in the model galaxies would 
make the model B — R colours agree with observations) . The 
model SMGs have significant overlap in their observed op- 
tical colours with BX/BM galaxies: 45% of model SMGs 
have colours consistent with BX galaxies, whilst another 
10% have colours consistent with BM galaxies. This mix 
of BX/BM is comparable to the observational constraints 
from C05. 



3.6 Restframe Optical Properties 

If we turn to longer wavelengths, we can better test the rest- 
frame optical properties of galaxies predicted by the models 
against the observations. In Fig. [S] we compare the predicted 
K-h anA magnitu des of model SMGs with the photometry 
from lSmail et al.l (|2004j)- 

Since the stellar population models 
do not include emission line components, we have corrected 
the observed A'-band photometry for galaxies between z — 
1.8-2.8 (where the strongest optical emission line seen in 
SMGs, Ha, falls in K) for the Ha line emission assuming 
the median Ha equi valent width of EWrest (Ha) = 75 ±25 A 
measured in SMGs l|Swinbank et al.ll2004l ). The correction 
to the if-band photometry is ^ 10%. In contrast to the 
optical bands, as Fig. [6] shows, the if-band magnitudes of 
observed SMGs are brighter at a given 850/im flux density 
than all but the brightest model SMGs. However, this plot 
hides the redshift evolution of the iC'-band magnitude and 
so a more useful comparison comes from the evolution of 
K-hand flux as a function of redshift. As Fig. [6] also shows, 
at z ~ 2 the median observed A'-band magnitude for SMGs 
is = 20.0 ± 0.3, which is approximately two magnitudes 
brighter than the model prediction of = 22.0 ± 0.9 (for 
both model SMGs and rSMGs). 



However, at z ~ 2, the observed /('-band samples the rest- 
frame V-band which is dominated by young stellar popu- 
lations. So this comparison is still sensitive to the precise 
geometry and degree of any dust obscuration - where the 
models may lack sufficient detail. A more robust test of 
the predicted photometry from model galaxies comes from 
a comparison in the rest-frame Tf-band where the effects 
of dust extinction are much reduced. The rest-frame near- 
infrared can also provide a better estimate of the mass of 
the stellar population. At z ~ 2, the rest-frame ii'-band is 
redshifted to ~ 6/im, and as such the Spitzer Space Tele- 
scope has provided unique insight into the stellar masses of 
far-infrared luminous galaxies. 

In particular, iBorvs et al.l (|2005l ) (see also 
[Alexander et al.l |2008| ) estimated rest-frame 2.2/im lu- 
minosities for a sample of spectroscopically confirmed 
SMGs in the GOODS-N region from broad-band photom- 
etry covering UBVRIzJK-\-3.(3/4.5/5.8/8nm which spans 
the rest-frame UV, optical and near-infrared at the redshifts 
of their sample, z = 0.6-2.9. They derive Mk = -26.8 ±0.4 
from their sample of ten galaxies at z > 1.5. They then 
estimate a light-to-mass ratio for the stella r population of 
LkJM = 3.2 using the Starburst99 model (|Leitherer et al.l 
119991 ) based on an assumed mean age o f 200 Myr and 
a MiUer-Scalo IMF (|Mifler fc Scald 1 19791 ) and so infer 
a typical stellar m ass of log(M*) = 2.512 5 ^ ^^0- 
iBorvs et all l|2005h estimate that the dust extinction in the 
rest-frame K-hand is only ~ 0.2mag for their galaxies, and 
so do not correct for it in their stellar mass estimates. In 
comparison, our model SMGs have median stellar masses 
M^, ~ 2.ll5;j5 ^ 10^° Mp), up to an o rder of magnitude 
lower than inferred bv lBorvs et al.l (|2005l ) (we note that the 
stellar masses quoted here from both observations and the 
models include the mass of living stars plus remnants). 

However, since deriving the stellar mass is sensitive 
to the mass-to-light ratio (which depends on the assumed 
IMF, as well as the age of the stellar population and the 
dust extinction), a more robust comparison is to simply 
compare the observed mid-infrared fiux densities. In Fig. [7] 
we compare the predicted 5.8/im fiux density with obser- 
vations compiled fr o m deep Spitzer surveys l|Egami et al.l 
l2004l : lHainhne|[2007l : iPope et al.|[2006h . For 2 = 2-3, obser- 
vations suggests a median 5.8/im flux density of 35 ± 8^tJy. 
In comparison, the model SMGs have 5.8/im flux density of 
3.6l2'4MJyi thus under-predicting the 5.8/im flux density by 
a factor ~ lOx. Thus the low rest-frame if-band fluxes in 
the models compared to observations seem to point to the 
same conclusion about the stellar masses in the model being 
too low. 

We note that the rest-frame stellar mass-to-light ratio 
for the model SMGs (including dust extinction) is Ijk / M = 
4.8 ± 2.9. This is not greatly dissimilar to the value used 
bv iBorvs et al.l (|2005l ). but given the different IMFs, stel- 
lar ages and dust extinctions assumed, this apparent consis- 
tency may be somewhat fortuitous. 

This effect of the predicted 5.8/im fluxes being lower 
than observed values is not only limited to the model SMGs. 
We also use the model to select Lyman-Break Galaxies 
which have also had extensive mid-infrared follow-up. The 
median 5.8/im flux density of 72 spectroscopically confirmed 
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Figure 6. Predicted observed X-band magnitude distribution as a function of redshift and 850/^m flux density for model SMGs and 
rSMGs compared to the observations (median values are plotted for the models). The error- bars map the central 67% of the distribution 
from GRASIL and are shown on alternate points for clarity. Both plots show that the predicted photometry of model SMGs under-estimates 
that measured from observations by ~ 2 magnitudes (~ 6x). 



LBGs between z = 1.5-3 from lShapley et al.1 (120051 ) is Ss.g = 
S^i's^Jy (including non detections). In contrast the median 
5.8/im flux density for model LBGs selected using the same 
colour cuts and redshift distribution is Ss.g ~ 0.4tQ'2/iJy, 
approximately a factor 4x lower than observations suggest. 

Are the low _ft'-band and 5.8/im flux densities for model 
SMGs an effect of the stellar population models rather than 
an effect of the stellar masses (i.e. is there some way that 
the models can produce more rest-frame K-hanA light for 
the same stellar mass)? 

One potential route to reconciling the models with ob- 
servations is to introduce a full treatment of the thermally- 
pulsating asymptotic giant branch (TP- AGE) phase into the 
stellar p opulations. Recen t stella r population synthesis mod- 
eling by iMaraston et al.1 l|2006l ) has shown that TP-AGB 
stars in post-starburst galaxies can contribute significantly 
to the rest-frame K-hanA luminosity on ^ 200-1000 Myr 
timescales. Indeed, the predicted stellar light-to-mass ratio 
at 250 Myr for a s ingle stellar populati on is L_ff/M~ 6- 
8 iMarastonl [l99i : IMaraston et al.1 120061 ). This acts to re- 
duce the observationally estimated stellar mass by a factor 
of up to 3 X if the luminosity weighted stellar populations 
in the SMGs are dominated by stars with ages of ^ 200- 
1000 Myr. However, since sub-mm galaxies are likely to be 
sustained bursts, the dominant stellar population is likely 
much younger, and the mass-to-light ratio evolves as a func- 
tion of the burst age. We note that the median age of a burst 
in a sub-mm galaxy in the models is 55I3Q'' Myr (where age 
denotes the time elapsed since the beginning of the burst). 
To test the likely effect on the GRASIL luminosities we take 
the existing star-formation historie s and construct t he com- 
posite stellar population using the iMarastonI l|l998l ) stellar 
population models. We use a fiat IMF {x = 0) in bursts 
and standard IMF in the quiescent mode, as in our fiducial 
model. The predicted if-band and 5.8/im fiux densities are 
not significantly effected: the observed K-hanA fiux density 
increases by ~ 40% at z = 2 whilst the 5.8^m flux density 
increases by only ~ 20% on average (Fig. [7|. This suggests 
a median 5.8^m fiux density for model SMGs of --^ 4/iJy for 



the Maraston stellar population models compared to 3/iJy 
for the Padova models at z — 2-3, still a factor 8-10 x lower 
than that inferred observationally. Thus, although the TP- 
AGB phase may have a small effect in the SED modeling, for 
extended bursts this phase is unlikely to account for the low 
if-band and 5.8/im luminosities predicted for model SMGs. 

Another possibility to reconcile the models would be 
to reduce the burst lifetimes such that the rest-frame near- 
infrared light becomes dominated by TP-AGB stars. For 
short bursts (i.e. < lOMyrs), the younger luminosity- 
weighted age for the stellar population in the SMGs results 
in a much hi gher light-to-mass rat io, with Lif/M = 20-50 
at 15-30 Myr (|Maraston et al.ll2006l '). However, these shorter 
bursts are effectively ruled out since the amount of energy 
required to heat the dust in the sub-mm phase remains fixed, 
shorter burst durations result in hotter dust temperatures 
and therefore the 850/im sub-mm counts are not matched. 
This leads us to conclude that the deficit of stellar mass is 
a real effect and can not be easily explained by either the 
introduction of the TP-AGB phase into the stellar popula- 
tion models or by dramatically reducing the burst lifetime. 
We return to this in §5. 

One other potential contributor to the low predicted 
mid-infrared fluxes is the effect of dust extinction in the 
model. Model SMGs at z ~ 2 typically have extinctions 
of 4 mag in the observed i?-band, 2.4 mag in the observed 
if-band, and 1.6 mag at observer-frame 5.8^m. This means 
that the unextincted 5.8/im flux would be 4x larger than 
the extincted value. If the 5.8/im dust extinction could be 
drastically reduced in the model SMGs, while still retain- 
ing a large enough extinction in the rest-frame UV to cause 
most of the UV luminosity to be reprocessed by dust, then 
this could remove some of the discrepancy between pre- 
dicted and observed 5.8/im fluxes for SMGs, although a sig- 
nificant offset, ^ 2, would remain. However, we note that 
there is support for these high extinctions in the K-hanA, at 
least for the emission-line gas. Indeed, studies of the Balmer 
decrement in SMGs suggest extinctions of order y4„=2.9±0.5 
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Figure 7. Predictions for the 5.8^m flux density of model SMGs as a function of 850/.tm flux density and redshift. In both panels we 
show the trend of the predicted photometry from GRASIL and the changes which occur when the star-formation histories assume the 
stellar population synthesis from iMarastonI lll998|) which i nclud e a full treatment of TP- AGB stars. The co mparison sample comprise 
SMGs with (secure) spectroscopic redshifts from iHainlind l|2007l ) ; lEgami et al ] Hooi) and IPope et al Both plots show that the 

predicted 5.8/im flux density for model SMGs is systematically lower than the observations (by a factor lOx), and we suggest that this 
deficit is due to the lower-than-expected stellar mass for model SMGs (M* ~ 10^'^ M© whilst observations suggest more massive systems, 
Nh ~ 10" Mq) 



jTakata et al.|[200^ ). We do not explore this possibility fur- 
ther here, but defer this to a future paper. 

Finally, we note that while AGN are ubiquitous in 
SMGs, th eir typical con tribution to the bolom etric emission 
is <, 10% ([Alexander et al. 2005 : Menendez-Delmestre et al] 
I2OO7I : [Alexander et al. 2()od ) suggesting that any contami- 
nation to the rest-frame 2.2^m flux from AGN activity is 
minimal a nd unlikely t o affec t the observed stellar masses. 
Indeed, as iBorvs et al.l (|2005l ) point out, the shape of the 
spectral energy distribution is well-described by a stellar 
model and remove from their analysis the small number of 
galaxies for which AGN potentially dominate the rest-frame 
near-infrared. 



4 MASSES AND EVOLUTION OF SMGS 

4.1 Kinematics and Dynamical Masses 

Deriving dynamical masses for high redshift galaxies relies 
on measuring rotation curves and/or line widths and sizes 
for the line emitting regions. The most reliable line width 
estimates for sub-mm galaxies come from reso lved dynam- 
ics t r aced through millimetric CO emission (|Greve et al.l 
I2OO5I : iTacconi et ailbogd 120081) and nebular emission (such 
as Ha: ISwinbank et al.ll2004 I2OO6I '). Converting these Une 



widths into dynamical masses is particularly problematic 
at high redshift since the sizes of the galaxies are poorly 
constrained, even with HST resolution. We therefore be- 
gin by comparing the line widths of the model and ob- 
served SMGs. Based on observations of approximately 30 
SMGs, GHc = 170 ±30kms-^ and ago = 200 ±45 km s'^ 
have been derived dSwinbank et al.ll2004 iGreve et aLll2005l : 
iTacconi et al.ll2006l I2OO8I ). For the model galaxies the veloc- 
ity dispersions and sizes are calculated assuming (i) conser- 
vation of angular momentum during the collapse of the dark 
matter halo and baryons; (ii) the size of a stellar spheroid 



remnant produced by mergers or disk instabilities is deter- 
mined by virial equilibrium and energy conservation (see 
ICole et al]|2000l for a detailed discussion). Disks are assumed 
to have an exponential surface density profile with half mass 
radius rdisc, whilst the spheroid follows an r^'^^ law in pro- 
jection with half mass radius (in 3D) Vbuige- The mass dis- 
tribution in the halo and the length scale of the disk and 
bulge are assumed to adjust adiabatically in response to each 
other: for the disk the total angular momentum is conserved 
whilst for the spheroid rVc(r) is conserved at rtuige- For each 
galaxy this results in a size and velocity dispersion for the 
spheroid and a size and rotational velocit y for the disk and 
have been tested against observations bv lCole et all (|2000D 
for local disks and I Almeida et al] l|2007h for local spheroids. 

For both model SMGs and rSMGs we find am = 
160±30kms~i which is in good agreement with the obser- 
vational constraints. In order to constrain dynamical masses 
we also need to include the sizes. However, observationally 
these are poorly constrained, and usually are taken to be 4- 
8kpc (which is approximately the size seen in i/ST observa- 
tions and in resolved spectroscopic imaging). Using simple 
dynamical models, observations suggest d ynamical masses 
of or der ~ 2-5 x lO" M© within 4-8 kpc (|Swinbank et all 
I2OO6I ). 

Since GALFORM predicts the size and rotational ve- 
locity (or velocity dispersion) of the bulge and disk, it is 
also possible to crudely compare the predicted dynamical 
masses. In Fig. |5] we show the variation of projected 1- 
dimensional velocity with spatial scale for sub-mm galax- 
ies which show multiple components in their spatially re- 
solved spectra, cith er from CO, IFU or longslit obser vations 
(Tecza et al. 2004: , Swinbank et all [2004 . 2006; Gro ve et al.1 
I2OO5I : ITacconi et alJbood I2OO8I ) and overlay the predicted 



trend of rotational velocity versus size for model SMGs. Al- 
though it is clearly difficult to compare the velocity offsets 
between merging components and the rotational velocity of 
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Figure 8. Velocity offset versus spatial offset for the multiple 
component s within sub-mm galaxies from IFU, longslit or CO ob- 
serva tions fTccza ct al."2004; 'Gre veet ai]|2005l : ISwinbank et al.l 
|2004 2006; Tacconi ct al. 2006, 2001). We plot the average ra- 
dius and velocity of the disk for SMGs and rSMGs from GAL- 
FORM showing that the predicted dynamical masses of SMGs are 
in reasonable agreement with those inferred observationally. We 
also show the median halo circular velocity and radius for model 
SMGs (210±30kms-i and 110±20kpc respectively) at the edge 
of the plot. 

disks, since the predictions for model SMGs lie within the 
scatter of the observations, this suggests that the model and 
observational dynamical masses are in reasonable agreement 
(moreover, in an ideal system the expected difference be- 
tween the velocity offset between merging components and 
the rotational velocity of a disk is a factor ~ 2x). However, 
this comparison would benefit from improvements in both 
the observations and theoretical predictions (e.g. laser guide 
star adaptive optics integral field spectroscopy would trace 
the line emitting regions on sub-kpc regions) coupled with a 
theoretical study of the distribution of detectable emission 
line gas within dusty, merging syste ms at high redshift (e .g. 
from a direct numerical simulation: lOkamoto et ai]|2005l ). 

4.2 Gas Masses 

At its most basic level, the star formation in SMGs is 
fueled from the reservoirs of cold H2 gas in these sys- 
tems. Observations of emission from CO rotational transi- 
tions are the primary means of tracing (dipole-less) molec- 
ular hydrogen, H2. The conversion between CO line lu- 
minosity and total cold gas mass is usually expressed via 
Mgas/l^'co = aM0 (Kkms"Vc^)"^- The standard Milky 
Way ratio is M(H2-f He)/L'eo = 4.6M0 (Kkms-Vc^)"^- 
However, in the extreme environment near the center of 
an ultra-luminous galaxy, where the CO emission originates 
from an inter-cloud medium that is essentially volume fill- 
ing, rather than from clouds bound by self-gravity, the CO 
luminosity traces the molecular mass as a whole. Exten- 
sive high resolution CO mapping of local [z < 0.1) ULIRGs 
has been used to trace their gas reservoir sizes and dynam- 
ics to derive an approximate ratio of M_tf,+jfe/Lr7o ~ 0. 8 
iDownes fc SolomonI [l993 : ISolomon fc Vanden Boual2005h . 
Adopting this value to convert the CO emission line lumi- 
nosities of high redshift sub-mm galaxies suggest s gas masses 
3.0 ± 1.6 X 10^° M© within the central 2 kpc (jCreve et all 



I2OO5I : iTacconi et ahl |2006| . I2OO8I ). Note that if a > 1 and 
assuming random orientations then the H2 gas mass would 
exceed the dynamical mass on average. 

In the model SMGs, the median gas mass Mgas = 
iAt\l X 10^° Mq. a direct comparison of the gas masses 
of observed SMGs with those in the models is difficult 
due to the uncertainty in the conversion between CO(3- 
2) and CO(l-O) line lumi nosity as well as th e conversion 
between Ij'co ^-^i^ Mgas (|Hainline et al.l [2OO6I ). Neverthe- 
less, we can at least test what the value of a would be 
given the predicted gas masses and observed line luminosi- 
ties. The median line luminosity of >5 m jy SMGs from 
iGreve et al"] (|2005l ) and lTacconi et al.l (|2006l ) (including non 
detections) is ^'co = 3 ± 1 x 10^^° Kkms~^ pc~^. Assum- 
ing the gas mass of model SMGs and rSMGs of M^as = 
iAtll X 10^° Mq and using a = Mgas/Lco would re- 
quire a — l.l^Qg( K kms^^ pc~^)~'^ to reproduce the ob- 
served Ij'qq , which is in good agreement with the value used 
to observationally infer gas masses from CO line luminosities 
in similarly active systems. 

4.3 Halo Masses 

The clustering of sub-mm galaxies encodes important in- 
formation on the underlying dark matter distribution (as 
well as the distribution of SMGs within and between the 
halos). A measurement of clustering therefore provides key 
constraints on models of galaxy formation. Alt hough clus- 
tering measurements for SMGs are tentative, iBlain et al.l 
(|2004bl ) showed the potential of using precise spectroscopic 
redshifts to detect clustering in modest-sized samples, de- 
riving a crude constraint on their clustering length of ro — 
9.8 ±3.0 Mpc (Fig.m. Using the the Millennium Simulation, 
Almeida et al. (2008 in prep.) predict that the correlation 
function of model SMGs from the same model used here is 
close to a power law of the form ^ = (r/ro)^ over more than 
three decades in separation, with 7 — —1.94 ± 0.05 and a 
correlation length Vo — 8. 8 ±0.3 Mpc, in excellent agreement 
with observations. 

In Fig. [9] we show the clustering lengths of halos of 
mass log(Mh/M0) = 11, 12, 13 measured from the Millen- 
nium Simulati on for the evolution o f the dark matter in the 
ACDM model (ISpringel et al.ll2005h . For a clustering length 
of 9.8 ± 3.0 Mpc the correspon ding halo mass is predicted 
to be Mh = S.ltl l X 10^^ M0 l|Gao et al.ll2"005l ). This com- 
pares to the predicted halo mass for the model SMGs of 
Mhaio = 3.6ti;5 X 10^^ M0. These resuhs confirm the high 
masses of the SMG host halos in the model. Moreover, the 
close similarity between the halo masses predicted from their 
clustering and the actual halo masses of the model SMGs im- 
plies that any "merger bias" ha s only a weak effect on their 
clustering (|Percival et al.ll2003h . 

4.4 Comparison to other galaxy populations 

The GALFORM model allows us to compare the SMGs and 
rSMGs with other high-redshift populations. In this sec- 
tion we briefly discuss the properties of radio unidentified 
SMGs, as well as optically faint radio galaxies (OFRGs) 
which have been proposed to have a similar redshift distri- 
bution and bolometric luminosities as radio-detected sub- 
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Figure 9 . Comoying correlation lengtii for SMGs from 
iBlain et al.l l|2004lj 'l compared to model SMGs, as well 
as in contrast to other populations of low- and high- 
redshift galax ies (Bahcall & Soncira 1983: McCarthy et al.ll200ll: 
Croom et a ^ [2 001; Adcl bcrgcr & Steidcl 2000; Will mer et al.l 



IQQ^ r iOverzier et al.ll2003h . The dashed lines show the expected 



correlation length of dark matter halos as a function of mass and 
redshift, as measured from the Millennium Simulation (Gao et al. 
2005). A clustering length of tq = 9.8 ± 3.0 Mpc corresponds to 



a halo mass of Mi, 



'^h = S.llig X iU- 1V10. 

that directly predicted for the model SMGs, suggesting that any 
merger bias is low in our model. 
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mm g alaxies but with comparatively higher du st tempera- 
tures IChapman et al.ll2004l : [Blain et aL l2004al ) as opposed 



to cooler temperatures (or higher redshifts) for the radio 
undetected SMGs. 

To investigate the properties of the model radio un- 
detected SMGs and their relation to the model SMGs, we 
select model galaxies with Sgso > 5mJy and Si. 4 < 30/iJy. 
These have a median redshift oi < z >= 2.4, ctz = 0.5, 
peaking at somewhat higher redshift than the model SMGs 
{<z>= 2.0). The bolometric luminosities of the model ra- 
dio undetected SMGs are comparable to the model SMGs 
(Lioi = LB^'j' * X 10^'^ L©), with characteristic dust temper- 
atures slightly colder than mo del SMGs (T^ = 28 ± 8K). 
This supports the suggestion of IChapman et al] (|2005l ) that 
the radio undetected SMGs are likely to represent the cold 
wing of the distribution of dust temperatures, but with a 
similar bolometric luminosity to the radio-detected SMGs. 
As noted before, this results in a modest bias against the 
highest redshift SMGs in radio-identified samples. 

The halo, stellar and gas masses of this population 
are Mhaio = 2.5+0^ x 10^^ M* = 1.5lo:9 x 10^° and 
Mgas = 3.41^;° X 10^° Mq respectively. Thus the properties 
of the model radio undetected SMGs significantly overlap 
those of model SMGs, with the main difference arising in 
their dust temperatures, where the model radio undetected 
SMGs have marginally colder dust temperatures and sim- 
ilar bolometric luminosities, resulting in a slightly higher 
redshift distribution compared to model rSMGs. 

Similarly, we investigate the properties of the model 
OFRGs: we define OFRGs as having Sgso < 3mJy, S1.4 > 
30AtJy and Iab > 22.5 to be consistent with lChapman et al.l 
l|2004h . We find that the model OFRGs have a much lower 



mean redshift than the the model SMGs, with < z >= 1.2, 
(Tz = 0.6. Since these galaxies are selected via their radio fiux 
density, their 850/^m/1.4GHz fiux density is lower than the 
model SMGs with S850/S1.4 — 17tg at 2; = 1 compared to 
model SMGs with Sg5o/Si.4 = 50^1° at the same redshift). 
The halo, stellar and gas masses of OFRGs are Stf x lO", 
4^3 x 10^ and l.5to,t ^ 10^° ^0 respectively, and bolomet- 
ric luminosities of order 6.5 ± 1.5 x 10^^ L©, approximately 
3-5 X lower than the SMGs. 

To relate the properties of model SMGs, radio unde- 
tected SMGs and OFRGs in Fig [10] we show the predicted 
850/im/1.4 GHz fiux density ratio as a function of 1.4GHz 
(rest-frame) radio luminosity, split into four redshift bins 
(z < 1.5, z — 2, z — 2.5 and z > 3). For comparison we 
overlay the SMGs from C05. We also show along the top axis 
the approximate bolometric dust luminosity corresponding 
to the radio luminosity, computed using the far-infrared- 
radio correlation defined by the model SMGs (which has 
the form log{Ltai/LQ) = -14.40 ± 0.15 + 1.12 ± O.O2L1.4). 

As expected, at z < 1.5-2 the model radio-undetected 
SMGs are cooler with lower bolometric luminosities (lower 
dust temperatures and lower bolometric luminosities com- 
bine to give the same apparent submm flux density as model 
SMGs). In contrast, at a > 2.5 the radio fiux density limit 
cuts into the typical dust temperatures of model SMGs, 
and so typical model SMGs (which would have been identi- 
fied at lower redshift) are no longer selected, leaving just 
the tail of model SMGs with the highest dust tempera- 
tures. The observations also show three very low redshift, 
cold sources (4% of the C05 sample) which are missed in 
the models. These might arise due to mis-identifications 
(e.g. 2 radio sources within the beam or an unrelated ra- 
dio source and a radio undetected SMG within the beam). 
It is also possible that the model is missing a component 
of exte nded cool ULIRG emission (so calle d cirrus galax- 
ies; e.g. lEfstathiou" fc Rowan- Robinsonlliooi ). We also note 
that there there are a few high redshift sources with low 
S'85o/'S'i.4 flux density ratios and high bolometric luminosi- 
ties, indicating possible contamination from AGN in some 
of the higher-redshift sources (or an additional radio com- 
ponent which is not included in the models). 

Fig. [To] also shows that the model OFRGs are predicted 
to have higher dust temperatures than model SMGs, and a 
range of bolometric luminosities from Lboi ~ 1O^^~^^'^L0. 

Finally, we also use the model to test what fraction 
of high-redshift ULIRGs are identified as SMGs. Selecting 
galaxies with bolometric luminosities comparable to SMGs 
(> 2 X 10^^ Lq), we find that the model SMGs comprise 
~ 50% and ~ 40% of the ULIRG population at 2 = 1.5-2.0 
and z = 2.5-3.5 respectively. 



4.5 Comparison between the z~4 and Z'^ 
populations 



2 SMG 



The negative K-correction experienced at 850/.im potentially 
allows very high redshi ft (up to g ^ 8) galax ies to be de- 
tected in the sub-mm (jBlain fc Longaij ll996D . The prop- 
erties of the first sub-mm galaxies (2 <; 4) and those at 
2 ~ 2 may be very different as the Universe doubles its age 
over this period. For this reason we have compared and con- 
trasted the predicted properties of the most distant sub-mm 
galaxies, at 2 ^ 4, with the more typical systems at 2 ~ 2. 
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Figure 10. Radio (1.4 GHz) luminosity versus far-infrared/radi o flux density ra t io as a function of redshift for model SMGs, radio 
undetected SMGs (ruSMGs) and OFRGs compared to SMGs from lChapman et al. For clarity we split the data into four redshift 

bins. We also show the approximate corresponding bolometric luminosity, computed using the far-infrared— radio correlation defined by 
the model SMGs. We also show the approximate characteristic temperature. The plot shows that at a fixed luminosity, the model radio 
undetected SMGs have systematically lower temperatures than the model SMGs at each redshift. We also note that (by selection) the 
model OFRGs have systematically higher temperatures at a fixed luminosity compared to the SMGs at all redshifts. 



In Fig[TT]we split the population into two redshift bins 
[z = 1-3 and z > 3.5) and compare the halo, stellar and gas 
mass functions. The model SMGs at z > 3.5 have typical 
halo, stellar and gas masses 2.5 x 10^^ M©, 1.0 x 10^° M© 
and 5x1O^''M0 respectively. Thus the the z = 1-3 counter- 
parts have halo and stellar masses which are a factor 1.25 
and 2x larger respectively, with the gas reservoirs a factor 
~2x smaller. The bolometric luminosities of the z — 1-3 
and z > 3.5 model SMGs are Lboi = l-8ti;o ^ lO^^L© and 
2.8li'3 xW^^Lq respectively. Thus the lower redshift model 
SMGs appear to have slightly lower (~ 1.5 x) bolometric lu- 
minosities, although there is significant overlap. This lumi- 
nosity evolu t ion su pports that claimed trend first noted by 
llvison et al.l l|2002l ) using deep radio and sub-mm imaging. 
We also note that at z = 4 the radio-detected fraction of 
SMGs is ^30%, which is much lower than at z = 2 where 
the radio detected fraction of SMGs is ~ 95%. 

Overall, these facts suggest that the z > 3.5 model 
SMGs are physically quite similar to model SMGs at z = 1- 
3. Since there is significant overlap in their properties, the 
rapid evolution in the redshift distribution and space den- 
sities of model SMGs between z ~ 2 and z 4 seems at 
first sight puzzling: the negative K-correction means that 
far-infrared luminous galaxies do not get fainter at 850/im 
as their redshift increases. The explanation of their rapid 
evolution lies in the hierarchical buildup of structure in the 
dark matter. The halo masses of model SMGs are similar at 
both epochs (> W^^Mq), yet the co-moving space density 




Figure 11. Cumulative mass functions for the halo, dynamical, 
stellar and gas masses for model SMGs split between at z ~2 
compared to those at z ~4. The figure shows that the z ^ 2 
model SMGs are more ubiquitous by a factor ~6x compared to 
those at z ~ 4, but their halo, dynamical, stellar and gas masses 
are only a factor ~ 20-40% larger suggesting that the SMG phase 
occurs in similar types of galaxy at both epochs. 



of > 10 Mq halos at z ~ 4 is five times less than at z 2. 
This evolution in the halo number density accounts for most 
of the change in the number density of model SMGs over this 
redshift range, although other factors also ent er, such as the 
fraction of halos undergoing major mergers (|Lacev fc Cold 
[l993), and the fraction of baryons contained in the form of 
cool gas. 
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5 DISCUSSION 

Our fiducial model has previously been shown to fit the prop- 
erties of local IRAS galaxies, as well as the number counts of 
faint sub-mm galaxies. Here, we have compared the SEDs, 
dynamical and halo masses of sub-mm galaxies against ob- 
servational data, showing good agreement with the values 
measured or inferred from observations. However, the largest 
discrepancy between the observations and the predictions 
from the model arises in the observed if -band and mid- 
infrared photometry. Although we find that in the model 
SMGs the recent starburst dominates the A"-band and mid- 
infrared fiux densities, contributing ~ 70% at 5.8/im and 
~ 80% in the A"-band, we find that the if-band fiuxes of 
model SMGs are under-predicted by a factor ~ 6x com- 
pared to observations, whilst the 5.8/im fluxes (which corre- 
sponds to rest-frame if-band at 2: ~ 2) are under-predicted 
by a factor ~ 10 x. This discrepancy is not limited to the 
model SMGs, but also applies to other model populations. 
For example, selecting Lyman-Break Galaxies (using ugr 
colours to select ubiquitous star-forming galaxies at 2: ~ 3) 
we find that the model LBGs at z = 2-3.5 have 5.8/im flux 
densities a factor ~ 4x lower than inferred from observa- 
tions. We show that including a full treatment of the TP- 
AGB phase in the stellar population modeling l|Marastonl 
1 19981 ) cannot make up for this deficit for any plausible star- 
formation history, since this increases the 5.8/im flux densi- 
ties by only ~ 20% on average. We also note that within the 
halo hosting the model SMG, 80-90% of the stellar mass is 
associated with the model SMG itself, thus including more 
of the stellar mass from satellites would have little effect on 
the integrated light. 

The low K-hsLTid and mid-infrared flux densities are 
a concern for this model, since this is a key indication of 
the stellar mass and dust reddening. The obvious inter- 
pretation is that the stellar masses of the model SMGs 
{Mstar ~2X 10^" Mq) are too low. However, in the model 
SMGs the total stellar mass formed in the current burst is 
^ 10% of the total stellar mass at the epoch of observa- 
tion (typically 1-4 x 10® Mq buih in a period of 55^30° Myr 
which is the median age of the burst in a model SMG) . Fur- 
thermore, the total stellar mass formed in all of the pre- 
vious bursts contributes less than 10% of the total stel- 
lar mass of the model SMG at the epoch of observation. 
This means that the burst mode is responsible for forming 
less than ^ 20% of the total stellar mass in model SMGs. 
Since the current burst dominates the luminosity, the rest- 
frame if-band light-to- mass ratio for the existing- and burst- 
stellar populations in the model SMGs are very different, 
with L/f/M=1.5to g and Lk/M=14.5±9.5 for the existing- 
and burst- stellar population respectively. Since the current 
burst is predicted to dominate the rest-frame K-hand lu- 
minosity, even though it contains only ~ 10% of the stellar 
mass one solutions would be to boost the mass in the cur- 
rent burst by a factor ~ 10. However, if everything else in the 
model was left unchanged, this would result in an increase 
in bolometric luminosity by the same factor. However, since 
the the bolometric dust luminosities of SMGs already agree 
with observations in the current model, this will violate the 
850/im counts. 

Our current model assumes a very flat IMF in bursts, 
with slope 2; = 0, so one could consider changing this slope. 



Would steepening the IMF slope to something closer to the 
Salpeter value x — 1.35 (so producing relatively more lower- 
mass stars) , while at the same time increasing the total burst 
mass, allow the model to increase the Tf-band luminosity 
while keeping the bolometric dust luminosity equal to ob- 
served values? To test this we examine the ratio of the K- 
band-to-bolometric luminosities for both a top-heavy and 
Salpeter IMF in bursts, but find that this only weakly de- 
pends on the IMF slope from a; = to a; = 1.5, presumably 
because both luminosities are dominated by massive stars 
at these young ages. Thus changing the IMF slope does not 
seem enough to boost the K-hand luminosity by a factor 
10 X for a fixed bolometric luminosity. 

One other potential avenue to increase the predicted 
near-infrared fiuxes in the model SMGs is to alter the dust 
extinction. As expected, the model SMGs have high dust 
extinctions, even in the rest-frame Tf-band, where the ex- 
tinction is typically about 1.6 mag. Therefore if extinction 
law could be altered in the model to greatly lower the extinc- 
tion in the rest-frame K-hand, while keeping the extinction 
high in the rest-frame UV, then the rest-frame Jf-band lu- 
minosities could be boosted by a factor ~ 4x while keeping 
the dust bolometric luminosities the same. This would go 
some way towards resolving the problem with the predicted 
5.8/im fiux densities being too low. However, we note that 
there is support for these high extinctions in the Tf-band, at 
least for the emission-line gas. Indeed, studies of the Balmer 
decrement in SMGs suggest y4„=2.9±0.5, thus removing all 
of the dust extinction in the observed Tf-band may violate 
observations constraints. We do not explore this possibility 
further here, but defer this to a future paper. 

Alternatively, is it possible to increase the stellar mass 
by condensing more gas into stars both prior to and during 
the sub-mm phase thus boosting the observed iC'-band and 
5.8/im luminosities? In order to prevent the over-cooling of 
baryons into stars (and therefore match the present day K- 
band luminosity function) the feedback prescription in the 
B05 model is particularly efficient, with an instantaneous 
feedback rate, /? = Mej/SFR= 1.1-2.0 (between z = 1-3) 
where M^j is the gas ejection rate. When recycling from 
dying stars is included this factor is even more extreme, es- 
pecially for a top-heavy IMF. In this case, the the recycled 
fraction is _R = 0.9, thus only 10% of the newly formed 
stellar mass ends up in long lived stars and remnants (thus 
the ratio of mass ejection rate to formation rate of long 
lived stars is 10-20 x during an SMG burst). This efficient 
feedback results in 55-65% of the baryonic mass in the hot 
phase for a model SMG at the epoch of observation, with 
a hot gas mass of Mhot ~ 8 x 10^" M© and a total baryonic 
mass of Mbary ~ 1-5 X 10"Mq. Thus ff the feedback could 
be reduced, and a substantial fraction of the mass contained 
within both the cold and hot gas reservoirs could be con- 
densed into stars either prior or during the burst, the result 
would be an increase in stellar mass of a factor ~ 6 x . This 
has the obvious draw-back that there would be no gas avail- 
able for future star-formation (unless the halo accretes gas 
from a merger). 

The most promising route to explaining both the near- 
and far-infrared luminosities in sub-mm galaxies may lie 
in combining presc riptions in the B05 model with those in 
iBower et all l|2006D which uses AGN feedback to quench gas 
cooling in high-mass halos and thus prevent the formation 
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of too many very high-mass galaxies by the present day, 
rather than relying on very strong supern ova feedback, as in 
the BOS model. In the lBower et alj (j2006l ) model, high-mass 
galaxies form earlier than in the BOS model, so the stellar 
meiss function at high redshift has more high-mass galax- 
ies (though both have very similar stellar mass functions 
at z = 0), and this gives better agreement with the ob- 
served evolution of the rest-fram e _R"-band luminosit y func- 
tion than for BOS. However, the iBower et al. model 
in its standard form with a single solar neighborhood IMF 
under-predicts the sub- mm cou nts by a factor 2 0x. We 
have investigated modifying the iBower et al. I l|2006l ) by m- 
troducing a top-heavy IMF in bursts, as in the BOS model. 
We find that we can obtain a comparably good fit to the 
8S0/im counts as the BOS model, but the same model then 
predicts a present-day K-hanA luminosity function which is 
too bright by about 1 magnitude at the bright end. How- 
ever, we have not explored the full parameter space of mod- 
els combining AGN feedback with a top-heavy IMF, and so 
it may be that a model can be found which explains both 
the number counts and rest-frame if-band luminosities of 
the SMGs, while also being consistent with the observed 
K-haxiA luminosity function at different redshifts. We will 
investigate this in more detail in a future paper. 



6 CONCLUSIONS 

In this paper we have used the GALFORM semi-analytic 
model for galaxy formation together with the GRASIL spec- 
trophotometric code to investigate in detail the prop erties 
of sub-mm galaxies in the model of lBaugh et al.ll|2005h . We 
mimic the observational selection in our model catalogue 
and identify galaxies both via their 8S0/im and 1.4 GHz 
(radio) emission. We find that the radio-identified SMGs 
(rSMGs) are predicted to make up approximately 75% of the 
whole population. This is in good agree ment with the frac- 
tion found in observational surveys (eg. Ilvison et al.ll2005l . 
[2003). The main difference between the radio-identified and 
radio-unidentified sub-mm galaxies is the characteristic dust 
temperature, where the radio undetected SMGs are ~ 10% 
cooler. This results in model SMGs having a median redshift 
of 2 = 2.0 whilst model rSMGs (which include a radio flux 
density threshold Si. 4 > 30/iJy) have a median redshift of 
z = 1.7. 

We show that the far-infrared colours of model rSMGs 
are consistent with observational data. Fitting the 350 and 
8S0/im photometry we show that the model far-infrared 
SEDs are well described by single component modified 
(/? = l.S) blackbodies with temperature of 32 ± SK (al- 
though GRASIL does not assume a single dust temperature 
either within or between galaxies). Integrating the model 
SED and the black-body fit, the inferred and true bolo- 
metric luminosities for the model SEDs also agree, with 
the single modified blackbody fit recovering an average of 
0.94±0.25 of the bolometric luminosity. However, the me- 
dian 850/im/1.4 GHz flux density ratio for SMGs predicted 
by the model is systematically lower than that observed by a 
factor ~ 1.26±0.24x , even after the model has been normal- 
ized to match the far-infrared-radio correlation for ULIRGs 
at 2 = 0. This could be due to evolution in the far-infrared- 



radio correlation, or due to modest contributions to the ob- 
served radio flux density from AGN. 

We show that the predicted velocity dispersions are in 
good agreements with observations, with am ~ 160 km s^^ 
for the model SMGs, compared to observational constraints 
of a = 170-200 km s~^ inferred from the kinematics of the 
redshifted Ha and CO emission. We also show that the pre- 
dicted gas masses are agree well with those derived from ob- 
servations using a — 0.8. Turning to the halo masses, we find 
that the model SMGs reside in halos with masses of M^aio ~ 
10^^ '*Mq, and a correlation length of = 8.8 ± 0.3 Mpc, 
in excellent agreement with th e tentative observati onal esti- 
mate of ro = 9.8±3.0 Mpc from lBlain et all l|2004bt ). We find 
no evidence for a significant merger bias in the clustering of 
SMGs in the model. 

However, we find that the predicted if-band and mid- 
infrared (5.8/im) flux densities of the model SMGs are up 
to a factor 10 x fainter than observed. We discuss the possi- 
ble reasons for this. The stellar masses of the model SMGs, 
M* ~ 10^° M0, are a factor 10 x smaller compared to ob- 
servational estimates based on the same mid-infrared data 
(M* ^ 10^^ M0), but the observational estimates are very 
sensitive to the assumed IMF, burst age and dust extinc- 
tion. The simplest explanation for this discrepancy in the 
near- and mid-infrared fluxes is that it is due to the stellar 
masses being too low in the model. However, other factors 
may be contributing, for example in the model SMGs, the 
dust extinction is quite large even at S.8^m, resulting in 
a 4x reduction in the observed flux. If the rest-frame visi- 
ble and near-infrared dust extinctions in the model SMGs 
were removed, this would improve (but not eliminate) this 
discrepancy (an offset of a factor 2x would remain). How- 
ever, we note that there is support for these high extinc- 
tions in the if-band, at least for the emission-line gas. We 
also discuss possible routes for increasing the near-infrared 
luminosities and stellar masse s, both within the context of 
the current (|Baugh et al.ll2005l ) model, which assumes a top- 
h eavy IMF in bursts, and using the prescriptions described 
in lBower et al.l ([2000), which uses AGN feedback to prevent 
too many baryons cooling i n massive halos but assumes a 
normal IMF throughout. The lBower et al.l (|2006l ) model pre- 
dicts more high-mass galaxies at high redshift, and a _Rr-band 
luminosity function in better agreement with observations 
at high redshift, but under-predicts the sub-mm counts by 
a factor 20. We will discuss the properties of s ub-mm galax- 
ies in mod els which combine fe atures from the lBaugh et al.l 
1 20051 ) and lBower et all (|2006h models in a future paper. 

We also inves tigate the evo l ution of the sub-mm galaxy 
population in the lBaugh et al.l (|200St ) model by comparing 
the halo and stellar masses for model SMGs at 2 = 1-3 to 
those at 2 > 3.5. The 2 > 3.S SMG population are predicted 
to make up approximately 2% of the total SMG population 
(with the 2 > 4 rSMGs contributing approximately 0.5% 
of the total population). Nevertheless, we find only a mild 
increase in the halo and stellar masses between 2^4 and 
2 ~ 2, and therefore attribute the strong evolution seen in 
the redshift distribution primarily to the evolution in the 
space density of massive halos hosting SMGs. 

As well as the need for the models to provide a bet- 
ter match to observations, a number of upcoming surveys 
and facilities will provide valuable observational constraints 
which can be used to further test them. In particular, deep 
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SCUBA2 surveys will probe both 450 and 850/im, providing 
much tighter constraints on the SEDs of galaxies with bolo- 
metric luminosities a factor 10 x lower than currently achiev- 
able. When combined with results from Herschel (which will 
probe 60-670/im), the temperatures, bolomotric luminosities 
and dust masses of SMGs will be constrained with much 
higher accuracy than currently possible from the sparse and 
low signal-to-noise 350 and 850pm photometry alone. Inter- 
estingly, the model predicts that there is a significant tail 
of faint 450pim selected galaxies at high rodshift (Laccy ct 
al. 2009 in prep). Selecting gala^cies with S450 > 2.5 mjy 
(the expected 5-cr flux density limit for SCUBA2) the model 
suggests that there should be > 8000, 850 and 150 galax- 
ies selected at 450/im per square degree at z > 2, 4 and 
5 respectively 50%, 5% and 1% of the total population 
respectively), suggesting that the fainter confusion limit at 
450/xm may provide the most efficient route to get an unbi- 
ased census for the far-infrared selected galaxies at the high- 
est redshifts, providing insights into the physical processes 
of star-formation occurring at the peak epoch of galaxy for- 
mation, « ~ 2. 



ACKNOWLEDGMENTS 

We would like to thank the anonymous referee for provid- 
ing a number of suggestions which significantly improved 
the content and clarity of this paper. We thank Claudia 
Maraston for providing us with stellar population models 
with a top-heavy IMF, and John Helly and Liang Gao for 
calculating the merger rates and clustering lengths of ha- 
los from the Millennium Simulation. We also gratefully ac- 
knowledge Dave Alexander, Richard Bower, Reinhard Gen- 
zel and Linda Tacconi for valuable discussions. AMS ac- 
knowledges an STFC fellowship. IRS & CMB acknowledge 
the Royal Society. This work was also supported by the 
STFC rolling grant for extragalactic astronomy and cosmol- 
ogy at Durham. 

REFERENCES 

Adelberger, K. L. & Steidel, C. C. 2000, ApJ, 544, 218 
Alexander, D. M., Bauer, F. E., Brandt, W. N. Horn- 

schemeier, A. E., Vignali, C, Garmire, G. P., Schneider, 

D. R, & Chartas et al. 2003, AJ, 125, 383 
Alexander, D. M., Brandt, W. N., Small, I., Swinbank, 

A. M., Bauer, F. E., Blain, A. W., Chapman, S. C, & 

Coppin et al. 2008, AJ, 135, 1968 
Alexander, D. M., Small, I., Bauer, F. E., Chapman, S. C, 

Blain, A. W., Brandt, W. N., & Ivison, R. J. 2005, Nature, 

434, 738 

Almeida, C, Baugh, C. M., & Laccy, C. G. 2007, MNRAS, 
376, 1711 

Bahcall, N. A. & Soneira, R. M. 1983, ApJ, 270, 20 
Barger, A. J., Cowie, L. L., Sanders, D. B., Fulton, E., 

Tariiguchi, Y., Sato, Y., Kawara, K., & Okuda, H. 1998, 

Nature, 394, 248 
Baugh, C. M. 2006, Reports of Progress in Physics, 69, 

3101 



Baugh, C. M., Lacey, C. G., Frenk, C. S., Granato, G. L., 

Silva, L., Brcssan, A., Benson, A. J., & Cole, S. 2005, 

MNRAS, 356, 1191 
Benson, A. J., Bower, R. G., Frenk, C. S., Lacey, C. G., 

Baugh, C. M., & Cole, S. 2003, ApJ, 599, 38 
Benson, A. J., Lacey, C. G., Baugh, C. M., Cole, S., & 

Frenk, C. S. 2002, MNRAS, 333, 156 
Blain, A. W., Barnard, V. E., & Chapman, S. C. 2003, 

MNRAS, 338, 733 
Blain, A. W., Chapman, S. C, Small, I., & Ivison, R. 2004a, 

ApJ, 611, 52 
— . 2004b, ApJ, 611, 725 

Blain, A. W. & Longair, M. S. 1996, MNRAS, 279, 847 

Blain, A. W., Small, I., Ivison, R. J., & Kncib, J.-P. 1999, 
MNRAS, 302, 632 

Borys, C, Small, I., Chapman, S. C, Blain, A. W., Alexan- 
der, D. M., & Ivison, R. J. 2005, ApJ, 635, 853 

Bower, R. G., Benson, A. J., Malbon, R., Hclly, J. C, 
Frenk, C. S., Baugh, C. M., Cole, S., & Laccy, C. G. 2006, 
MNRAS, 659 

Bressan, A., Silva, L., & Granato, G. L. 2002, A&AP, 392, 
377 

CariUi, C. L. & Yun, M. S. 1999, ApJL, 513, L13 
Chapman, S. C, Blain, A. W., Ivison, R. J., & Small, I. R. 

2003, Nature, 422, 695 

Chapman, S. C, Blain, A. W., Smail, I., & Ivison, R. J. 

2005, ApJ, 622, 772 

Chapman, S. C, Smail, I., Blain, A. W., & Ivison, R. J. 

2004, ApJ, 614, 671 

Clements, D. L., Vaccari, M., Babbedge, T., Oliver, S., 
Rowan- Robinson, M., Davoodi, P., & Ivison et al. 2008, 
MNRAS, 387, 247 

Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, 
MNRAS, 319, 168 

Condon, J. J. 1992, ARAA, 30, 575 

Coppin, K., Chapin, E. L., Mortier, A. M. J., Scott, S. E., 
Borys, C, Dunlop, J. S., Halpern, M., & Hughes ct al. 

2006, MNRAS, 372, 1621 

Coppin, K., Halpern, M., Scott, D., Borys, C, Dunlop, J., 
Dunne, L., Ivison, R., & Wagg et al. 2008, MNRAS, 123 

Coppin, K. E. K., Swinbank, A. M., Neri, R., Cox, P., 
Smail, I., Ellis, R. S., Geach, J. E., & Siana et al. 2007, 
ApJ, 665, 936 

Cowie, L. L., Barger, A. J., & Kneib, J.-P. 2002, AJ, 123, 
2197 

Croom, S. M., Shanks, T., Boyle, B. J., Smith, R. J., Miller, 

L., Loaring, N. S., & Hoyle, F. 2001, MNRAS, 325, 483 
Downes, D. & Solomon, P. M. 1998, ApJ, 507, 615 
Dunne, L., Eales, S. A., & Edmunds, M. G. 2003, MNRAS, 

341, 589 

Bales, S., Bcrtoldi, F., Ivison, R., CariUi, C, Dunne, L., & 

Owen, F. 2003, MNRAS, 344, 169 
Efstathiou, A. & Rowan-Robinson, M. 2003, MNRAS, 343, 

322 

Egami, E., Dole, H., Huang, J.-S., Perez-Gonzalez, P., Le 
Floc'h, E., Papovich, C, & Barmby, P., et al.,. 2004, 
APJS, 154, 130 

Gao, L., Springel, V., & White, S. D. M. 2005, MNRAS, 

363, L66 

Genzel, R., Baker, A. J., Tacconi, L. J., Lutz, D., Cox, P., 

Guilloteau, S., & Omont, A. 2003, ApJ, 584, 633 
Granato, G. L., Lacey, C. G., Silva, L., Bressan, A., Baugh, 



The Properties of Sub-mm Galaxies in Hierarchical Models 17 



C. M., Cole, S., & Frenk, C. S. 2000, ApJ, 542, 710 
Greve, T. R., Bertoldi, F., Small, I., Nerl, R., Chapman, 

S. C, Blain, A. W., Ivison, R. J., & Genzel et al. 2005, 

MNRAS, 359, 1165 
Greve, T. R., Ivison, R. J., Bertoldi, F., Stevens, J. A., 

Dunlop, J. S., Lutz, D., & Carilll, C. L. 2004, MNRAS, 

354, 779 

Greve, T. R., Pope, A., Scott, D., Ivison, R. J., Borys, C, 
Conselice, C. J., & Bertoldi, F. 2008, ArXiv e-prints, 806 

Mainline, L. J. 2007, Ph.D. Thesis, Caltcch Astronomy 

Hainline, L. J., Blain, A. W., Greve, T. R., Chapman, S. C, 
Smail, I., & Ivison, R. J. 2006, ApJ, 650, 614 

Helou, G., Soifer, B. T., & Rowan-Robinson, M. 1985, 
ApJL, 298, L7 

Hughes, D. H., Serjeant, S., Dunlop, J., Rowan- Robinson, 

M., Blain, A., Mann, R. G., Ivison, R., & Peacock et al. 

1998, Nature, 394, 241 
Ivison, R. J., Grevo, T. R., Dunlop, J. S., Peacock, J. A., 

Egami, E., Smail, I., & Ibar et al. 2007, MNRAS, 380, 199 
Ivison, R. J., Grove, T. R., Smail, I., Dunlop, J. S., Roche, 

N. D., Scott, S. E., Page, M. J., & Stevens et al. 2002, 

MNRAS, 337, 1 
Ivison, R. J., Smail, I., Bargcr, A. J., Kncib, J. -P., Blain, 

A. W., Owen, F. N., Kerr, T. H., & Cowie, L. L. 2000, 

MNRAS, 315, 209 
Ivison, R. J., Smail, I., Dunlop, J. S., Greve, T. R., Swin- 

bank, A. M., Stevens, J. A., Mortier, A. M. J., & Serjeant 

et al. 2005, MNRAS, 364, 1025 
Ivison, R. J., Smail, I., Le Borgne, J.-F., Blain, A. W., 

Kneib, J.-P., Bezecourt, J., Kerr, T. H., & Davies, J. K. 

1998, MNRAS, 298, 583 
Kennicutt, R. C. 1998, ARAA, 36, 189 
Kovacs, A., Chapman, S. C, Dowell, C. D., Blain, A. W., 

Ivison, R. J., Smail, I., & Phillips, T. G. 2006, ApJ, 650, 

592 

Lacey, C. & Cole, S. 1993, MNRAS, 262, 627 

Lacey, C. G., Baugh, C. M., Frenk, C. S., Silva, L., 
Granato, G. L., & Bressan, A. 2008, MNRAS, 385, 1155 

Leitherer, C, Schaoror, D., Goldadcr, J. D., Delgado, 
R. M. G., Robert, C, Kunc, D. F., dc Mello, D. F., De- 
vost, D., & Hockman, T. M. 1999, APJS, 123, 3 

Lilly, S. J., Eales, S. A., Gear, W. K. P., Hammer, F., Le 
Fevre, O., Crampton, D., Bond, J. R., & Dunne, L. 1999, 
ApJ, 518, 641 

Maraston, C. 1998, MNRAS, 300, 872 

Maraston, C, Daddi, E., Renzini, A., Cimatti, A., Dickin- 
son, M., Papovich, C, Pasquah, A., & Pirzkal, N. 2006, 
ApJ, 652, 85 

McCarthy, P. J., Carlbcrg, R. G., Chen, H.-W., Marzke, 
R. O., Firth, A. E., Elhs, R. S., Persson, S. E., & McMa- 
hon et al. 2001, ApJL, 560, L131 

Menendez-Delmestre, K., Blain, A. W., Alexander, D. M., 
Smail, I., Armus, L., Chapman, S. C, Prayer, D. T., Ivi- 
son, R. J., & Teplitz, H. I. 2007, ApJL, 655, L65 

Miller, G. E. & Scalo, J. M. 1979, APJS, 41, 513 

Okamoto, T., Eke, V. R., Frenk, C. S., & Jenkins, A. 2005, 
MNRAS, 363, 1299 

Overzier, R. A., Rottgering, H. J. A., Rengelink, R. B., & 
Wilman, R. J. 2003, A&AP, 405, 53 

Percival, W. J., Scott, D., Peacock, J. A., & Dunlop, J. S. 
2003, MNRAS, 338, L31 

Pope, A., Scott, D., Dickinson, M., Chary, R.-R., Morrison, 



G., Borys, C, Sajina, A., & Alexander, D. et al. 2006, 
MNRAS, 370, 1185 
Salpeter, E. E. 1955, ApJ, 121, 161 

Sanders, D. B., Mazzarella, J. M., Kim, D.-C, Surace, 

J. A., & Soifer, B. T. 2003, AJ, 126, 1607 
Scott, S. E., Fox, M. J., Dunlop, J. S., Serjeant, S., Peacock, 

J. A., Ivison, R. J., Oliver, S., & Mann, R. et al. 2002, 

MNRAS, 331, 817 
Shapley, A. E., Steidel, C. C, Erb, D. K., Reddy, N. A., 

Adelberger, K. L., Pettini, M., Barmby, P., & Huang, J. 

2005, ApJ, 626, 698 

Silva, L., Granato, G. L., Bressan, A., & Danese, L. 1998, 
ApJ, 509, 103 

Smail, I., Chapman, S. C, Blain, A. W., & Ivison, R. J. 

2004, ApJ, 616, 71 
Smail, I., Ivison, R. J., & Blain, A. W. 1997, ApJL, 490, 

L5 

Smail, I., Ivison, R. J., Blain, A. W., & Kneib, J.-P. 2002, 

MNRAS, 331, 495 
Smail, I., Ivison, R. J., Owen, F. N., Blain, A. W., & Kneib, 

J.-P. 2000, ApJ, 528, 612 
Solomon, P. M. & Vanden Bout, P. A. 2005, ARAA, 43, 

677 

Springel, V., White, S. D. M., Jenkins, A., Frenk, C. S., 
Yoshida, N., Gao, L., Navarro, J., & Thacker et al. 2005, 
Nature, 435, 629 

Swinbank, A. M., Chapman, S. C, Smail, I., Lindner, C, 
Borys, C, Blain, A. W., Ivison, R. J., & Lewis, G. F. 

2006, MNRAS, 371, 465 

Swinbank, A. M., Smail, I., Bower, R. G., Borys, C, Chap- 
man, S. C, Blain, A. W., Ivison, R. J., & Howat, S. R. et 
al. 2005, MNRAS, 359, 401 

Swinbank, A. M., Smail, I., Chapman, S. C, Blain, A. W., 
Ivison, R. J., & Keel, W. C. 2004, ApJ, 617, 64 

Tacconi, L. J., Genzel, R., Smail, I., Neri, R., Chapman, 
S. C, Ivison, R. J., Blain, A., & Cox et al. 2008, ArXiv 
e-prints, 801 

Tacconi, L. J., Neri, R., Chapman, S. C, Genzel, R., Smail, 
I., Ivison, R. J., Bertoldi, F., & Blain et al. 2006, ApJ, 640, 
228 

Takata, T., Sokiguchi, K., Smail, I., Chapman, S. C, 

Geach, J. E., Swinbank, A. M., Blain, A., & Ivison, R. J. 

2006, ApJ, 651, 713 
Tecza, M., Baker, A. J., Davies, R. I., Genzel, R., Lehnert, 

M. D., Eisenhauer, F., Lutz, D., & Nesvadba et al. 2004, 

ApJL, 605, L109 
Vlahakis, C, Eales, S., & Dumic, L. 2007, MNRAS, 379, 

1042 

Webb, T. M. A., Lilly, S. J., Clements, D. L., Eales, S., 
Yun, M., Brodwin, M., Dunne, L., & Gear, W. K. 2003, 
ApJ, 597, 680 

Willmer, C. N. A., da Costa, L. N., & Pellegrini, P. S. 1998, 
AJ, 115, 869 

Younger, J. D., Fazio, G. G., Huang, J.-S., Yun, M. S., 
Wilson, G. W., Ashby, M. L. N., Gurwell, M. A., & Lai 
et al. 2007, ApJ, 671, 1531 

Yun, M. S., Reddy, N. A., & Condon, J. J. 2001, ApJ, 554, 
803 

This paper has been typeset from a T^/ file prepared 

by the author. 



